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' This is an analysis of the itinerant model for antiferromagnetism, in which is included both on-site 

and inter-site electron correlations. We also consider the band degeneration, which brings into the 
■ Hamiltonian the on-site exchange interactions. 

' The Green function technique is used and the coherent potential approximation (CPA) decoupling 

for the on-site Coulomb repulsion. This decoupling combined with the modified Hartree-Fock ap- 
proximation for the inter-site interactions generates the change in shape of the spin bands with 
growing interaction constants, which is described by the correlation factors and which decreases the 
kinetic energy of the system. The effective Hartree field and the gain in kinetic energy due to the 
on-site and inter-site correlation factors drive the antiferromagnetism. The on-site and inter-site in- 
^ , teractions act together towards the antiferromagnetism. Their cooperation decreases the interaction 

?— j ■ constants required for the antiferromagnetic ordering. 

This new approach allows for the antiferromagnetic instability in the purely itinerant model at the 
half- filling in the split band limit. This situation describes the high temperature superconducting 
cuprates. 

s 

o 

1. Introduction 

(N 

^ ■ The basic model for magnetic ordering of itinerant electrons in solids is the Hubbard model [1]. The largest 
interaction between itinerant electrons is the on-site Coulomb repulsion, U = (i, i | 1/r \ i, i) , where i is the lattice 
site index. In the mean-field approximation the Hubbard model leads to the well-known Stoner model for magnetism. 

Extension of the Hubbard model is the model in which in addition to the on-site Coulomb repulsion U, there is 
also the Coulomb exchange on-site (between different orbitals) interaction «/„. The magnetic solutions of the three- 
dimensional Hubbard model did not make the exchange field necessary for antiferromagnetism to be small enough to 
be justified physically (see e.g. Ref. [2]), therefore the authors failed to arrive at the magnetic state. More recently 
Hirsch [3,4] and others (see e.g. Refs [5-7]) also included the nearest-neighbor inter-site interactions; | 1/r | k,l), 
where (k,l) = and i,j are the nearest-neighbor lattice sites. 

In our degenerate band model, we introduced into the Hamiltonian the on-site Coulomb repulsion 
. U = (iX, iX | l/r\i\,i\), the on-site two sub-band interactions (iX,iv\l/r\i'd,ivj) (where the sub-band indices 
(zu,"3) = (X,f) and A ^ v), as well as the inter-site nearest-neighbor interactions | 1/r \ k,l) (for (k, I) = 
O ■ and i ^ j). 

Through analyzing the inter-site correlations, the product of four operators is dealt with; two of which are replaced 
£> . by its stochastic value. The stochastic value of the two-operator product has been interpreted as the product of 
probabilities for electron transfer between two atoms. This interpretation explains that Hirsch's [4] the average bond 
occupation for spin a; i„, comes from the standard probabilities used in the coherent potential approximation (CPA) 
method [8] when it is applied to the inter-site correlation in the lowest mean-field (or Hartree-Fock) approximation. 
This connection will allow the expansion of the model in the future to more realistic cases of inter-site interactions of 
arbitrary strength. In this paper the weak inter-site interaction is accompanied by the weak or strong on-site Coulomb 
interaction. In the latter case the self-energy has higher order expansion in powers of interaction constant over the 
bandwidth. The spin band narrowing (broadening) is arising from inter-site interactions (see Refs [3,9]). 

The paper is organized as follows. In Section 2, the model Hamiltonian is put forward and in Section 3 the 
CPA formalism is developed to treat the inter-site Coulomb correlation. In Section 4, we set up the model for 
antiferromagnetism, which includes the Coulomb on-site and inter-site correlations. Numerical examples are presented 
in Section 5 based on the constant density of state (DOS) for the weak inter-site interactions and weak or strong 
on-site Coulomb interaction. Finally, Section 6 is devoted to the conclusions. 
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2. The model 



The Hamiltonian for the single degenerate band can be written in the form which has been introduced previously 
(see Ref. [10]) 

H = -t ^2 ( c ia c 3° + h - c ) - Vo^ht - i^°^<n JCT + U^2n ia ni- a 

<ij><7 i i.er ia 

+ J E C t C jo' C ™' C j<r +V E «*»J + J ' E C t C i-* C j-? C 3<T 
<ij>aa' <ij> <ij>a 

+ ]T At(cfc ja + h.c.)(h i - <r + n j - <r ) (1) 

<ij>cr 

where to is the nearest-neighbors hopping integral, fi is the chemical potential, cf a (c i(7 ) creates (destroys) an electron 
of spin a on the ith lattice site, fii a = c\ a Ci a is the electron number operator for electrons with spin a on the ith lattice 
site, hi — hia + hi-u is the operator of the total number of electrons on the ith lattice site, nf is the average number 
of electrons on sites i with spin a. The model Hamiltonian (1) includes explicitly the band degeneration already in 
the Hartree-Fock approximation. The internal field constant F® D is given by 

F? n = (p-l)(J in + j! n + V in ), (2) 

where p is the number of sub-bands in the band and J m = (iX,iv \ 1/r \ iv,iX), J/ n = (iX,iX \ 1/r \ iv,iv), 
Vin = (iX, iv | 1/r | iX, iv) are the on-site two sub-band inter orbital interactions of the following types: exchange, pair 
hopping and density-density, respectively (see Ref. [10]). The indices A, v numerate the sub-bands in the degenerate 
single band. 

The remaining potential part of the Hamiltonian contains the on-site Coulomb repulsion U = (iX, iX | 1/r | iX, iX), 
the inter-site exchange interaction J = (i, j \ 1/r \ density-density interaction V = \ 1/r | i, j), pair-hopping 
interaction J' = | 1/r | j, j) and assisted hopping interaction At = \ 1/r | j, i). 

The on-site Coulomb repulsion will be set aside for the time being and only inter-site terms will be consider in the 
Hamiltonian (1). The following form is arrived at 

H = ^2(e° k - Ho)h ka - F? n ^ n ifiia + Y A *("*-<7 + nj-v) ( C ia C ^ + C t<r Ci °) 

ka ia <ij>a 

+ E [ Jc t-^- + J'4-aC 1 -<r + (J-V)c+ a C la ]c+C la , (3) 

<ij>0 

where h ka = c ka Ck<j, c\ a {c ka ) creates (destroys) an electron of spin a and momentum k and e° is the unperturbed 
dispersion relation for electron 

e° k = -ton, n= E (4) 

<i,j> 

Next we will analyse Hamiltonian (3) using the coherent potential (CP) approximation scheme and its lowest order 
mean field decoupling. 



3. CP approximation for inter-site interactions 

The CPA idea is now applied to the inter-site weak interactions; J, J', V and At, in Hamiltonian (3) and the 
stochastic value will replace each operator product in the square bracket. For example, the operator product; c^ a Ci a 

will be replaced by its stochastic value c^ a c ia . It takes the values 1 or 0. The probability of value 1 is the probability 
of electron with spin a hopping from the i to the j lattice site and back. This probability denoted by l a below, is 
given by the average of two products. One is the product of probabilities that there is an electron with spin a on the 
i site and that the j site has empty states; nf (n* CT — nj)/ (nf + (n 1 - 7 — rrj )), and the second one is the probability 
of the opposite jump; nj(n- cr — nf)/ (rrj + (n t i a — nf)). The quantity nf is the average number of electrons on sites 
i with spin a and n* CT is the total capacity of the sub-band a on site i. The total hopping probability, called is 
given by the average of above probabilities 



2 



nf(nf -nj) 



For the weak correlation, when the band is not split by the Coulomb repulsion, we have n* CT = 1 and 



(5) 



I. = 



This expression for the ferromagnctism (nf 



<(!-<) 



nf(l-nf) 



< + (l-nJ) nj + (l-nf) 

= rij =n cr ) can be simplified to 1 
% a =<(!-<). 



(6) 



(7) 



In the case of strong correlation when the band is split into the lower and upper Hubbard sub-band, n* CT is the total 
capacity of the sub-band in which the chemical potential is located. For example, for the lower Hubbard sub-band 
we have; n- CT = 1 — n~ a . These values of n' CT together with the definitions of nf for antifcrromagnetism will be used 
later to calculate I a = I- a = Iaf in the case of antiferromagnetism. 

After adding up and subtracting the inter-site self-energy S \ 2 into the Hamiltonian (3) and performing the above- 
described stochastic approximation, this Hamiltonian takes on the following form: 



H = E - + E ( £ ~ - {c+c ia + c+ c i(7 ) + {M? - F&n<{) n lff , 

kcr (ij) ia 

cr 

where the effective dispersion relation is given by 

£ k — e k + 7fc S l,2 = £ k L o- 

The quantity L a is the parameter of the bandwidth change taking the following form: 

^1,2 



L D ° 1 



to 



(8) 



(9) 



(10) 



where Dq = zto is the unperturbed half bandwidth, z is the number of nearest neighbors and D a is the spin- 
dependent effective (perturbed) half bandwidth, for which we can write D a — zt a , where the perturbed hopping 
integral t a — L a t$. 

The quantity Mf in Eq. (8) is the stochastic energy describing the inter-site interactions 



(11) 



where is the sum over the nearest neighbors of the lattice site i, the bar over the operators product defines this 
product as the stochastic variable. 

The stochastic energy M? together with the on-site Coulomb repulsion can be subsequently analyzed in the CP 
approximation (see Ref. [8]). 

The stochastic potential in Hamiltonian (8) can be expressed as 



£ = ( J + j')c+_ a Ci- a + ( J - V)c+ a C la + At (m- a + Uj-a) , 



(12) 



with the stochastic values and 1 for each operator product in Eq. (12). The stochastic potential e will take on 16 
different values, ej, with corresponding probabilities, Pf , given by Eq. (A2) in Appendix. For these probabilities and 
potentials it is now possible to write the following equation for the self-energy SJ 2 , which describes only the inter-site 
interactions in Hamiltonian (8) 



1 Micnas et al. [9] and Hirsch [4] define the average band occupation for spin a as i CT = j jf ( — 75^) f r7 (e)de. For the constant 

-Do 

DOS this will give i a = nf (1 — nj) which is identical to our l a from Eq. (7). 
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where the Slater-Koster function (e) has the following form 
The Slater-Koster function can be written as 



N , 



where the Green function G a (e, k) is given by 



G a (e, fc) = „— = -. (16) 

In the first order approximation, for the weak inter-site interaction, we obtain the following formula for the inter-site 
self-energy EJ 2 

16 

^1,2 = E £ i p i = ( J - w + ( J + J ') J - + At ( n r ff + «7 CT ) (17) 

;=i 

and the parameter of the bandwidth change in Eq. (16) is given by 



l. . i - E* _ , _ ^<„ r + „ r) _ _ (iz^, (18) 

to to to to 



where nf is the average number of electrons with spin a on the zth lattice site. 

In the case of ferromagnetism, the parameter of the bandwidth change L a was spin dependent (see e.g. Refs [3,10]) 
which caused the spin dependence of the hopping integral t a . The value of the hopping integral for the majority spin 
a decreased with respect to the hopping integral for spin —a. 

For the antiferromagnetism the indices i,j belong to the neighboring sub-lattices a, (3 with opposite magnetic 
moments. In the result parameter; I a = I- a = Iaf is spin independent (i.e. it does not depend on the first power of 
antifcrromagnetic moment m) and the bandwidth factor from Eq. (18) 

L AF = ^- = l-n^ + ^(-2J-J' + V)I AF (19) 

L>0 to to 

is also spin independent (here n is the total number of electrons). 

The factor Laf is spin independent and because of it the dispersion relation will also not depend on the spin 
direction e£ = £,7°" = Sk = £°£af- The Green function from Eq. (16) will now take on the form 

G{e',k) = - r ^—, e' = £ + ^ , (20) 

which will be used in the next section in the case of antifcrromagnetic ordering. 

In the first order approximation the stochastic inter-site energy Mf given by Eq. (11) becomes the molecular field, 
which is given by the expression: 

M? = -J^2'n° + zVn + zAtI AF . (21) 

j 

In further analysis of the antiferromagnetism, the first order approximation for the parameter of the bandwidth 
change, Laf will be used (Eq. (19)) and the stochastic energy M? will be given by Eq. (21). 
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4. Antiferromagnetism (AF) 



For the antiferromagnetism, the diagonalization of Plischke and Mattis [11], Brouers [12] and Mizia [13] will be 
used. The crystal lattice will be divided into two interpenetrating sub-lattices a, (3, with the average electron numbers 
equal to 

4° = ^, (22) 

where m is the antiferromagnetic moment per atom in Bohr's magnetons. 

The reduced part of the Hamiltonian (8) together with the on-site Coulomb correlation is given by 

H=-tJ2(a+(3 ja +(3+a i(T ) - M o E n? 

i.j «-7,c 

+ E (M^-F° n n;)nl + U ]T n^_ ff (23) 



where a^ T (ai <J ) and (3f a {Pia) are the creation (annihilation) operators for an electron of spin a on sub-lattice a and 
(3 respectively, nf a — af^a^ is the electron number operator for electrons with spin a on the sub-lattice a, and 
t = toLxF is the effective hopping integral. The exchange potential in the mean-field approximation is different than 
in the case of ferromagnetism since the moments on neighboring lattice sites are opposite to each other and in effect 
both the inter-site and on-site component of the exchange potential will depend on the type of lattice site (a or (3). 
Combining Eq. (21) with Eq. (22) we obtain for the inter-site component that 

M? a{f)) = M° {0) = -Jzn% (a) + zVn + zAtI AF . (24) 

Now we will derive the equations for dispersion relation, particle number and magnetization for the antiferromagnetic 
state using Hamiltonian (23) and the Green function technique. Mean-field approximation is used for the parameter of 
the bandwidth change Laf and for the stochastic energy Mf = (Eq. (24)). This approximation was performed 

directly on the Hamiltonian, as there is no difference with the case when it is done within the Green function chain 
of equations. For the on-site interaction, U, we use the standard CPA, which for simplicity, in further analysis, is 
treated in the weak and strong correlation limit. For the weak correlation, (U <C D), we use first order approximation 
in interaction constant over the bandwidth. This is equivalent to the Hartree-Fock or virtual crystal approximation. 
The second case is the high correlation, (U ^> D), approximation. This approach will allow for easy extension of the 
calculations to the arbitrary strength of the on-site interaction later on. 

The main idea of the CPA formalism [8] is used now. We split the above stochastic Hamiltonian (23) into a 
homogeneous part 

i.j *,7,c icr i& 

(■y = a,/3) 

and a stochastic part 

Hi = E ^ + J2 (Vp - E2)nf„ (26) 



where 



H = Ho- zVn - 2zAtI AF (27) 



is the effective chemical potential, are the self energies on sites 7 = a((3) for electrons with spin a and V^iP) arc 
the stochastic potentials given by 

y<j = \ V la(f3) = - F iLK(0) - zJn l( a ) 

1 \f<? —TT—F°n' y — 7 Tn a ' 
[ V 2a(P) - U P in n a(p) ZJn f3(a) 

with probabilities 
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per _ -<J 



a(0) 



(28) 



The self-energies fulfill the equations 

E P n T^ 2 V = ' ^ = a ^ 

7 l-(V£-E-)*3» 

with the Slater-Koster function FJ 7 (e) in the following form 



(29) 



(30) 



and GJ y (e,k) given by Eq. (35) below. We transform Hamiltonian (25) into momentum space, and use it in the 
equations of motion for the Green functions: 



e((A;B)) E = ([A,B] + ) + (([A,H }_;B)) £ 
where (A, B) £ { a kai akc " @ka> @ k °)- ^ ne following equations are obtained 



e + fi - K -£fc 
-£ fc £ + /i - 



Gr(e,fc) G^(e,fc) 
G^(e,fe) G^(e,fc) 





'10" 




1 



(31) 



(32) 



with e.g. G™(e,k) = «a fctr ;a+ )), G*"(e,fc) - ((/3 fctr ; a+ )). 

Solving the set of these equations we arrive at the following expressions for the Green functions: 



G™(e,k) Gf{e,k) 
G^(e,k) G?/(e,k) 



1 

Dot 



e + n - e k 



where 



Dot = (e + ft - ££) (e + fi-'Sl)- e£. 



From Eq. (33) we have for 7 = a(/3) 



G^( £ ,fc) = - 



1 / e + ^- S flC„l 

^[G( £eff ,fc)-G(- £eff ,fc)], 



2\/ £ + /z-£ 



a(/3) 



with 



£eff = J(e + M - SS) (e + /x - S^) 



and G(£ e ff, &) given by Eq. (20) with e' replaced by e e g and /xo by it from Eq. (27): 

1 



G (e eff , k) 



e e ff — £fc 



e fe = e° k L AF . 



(33) 



(34) 



(35) 



(36) 



(37) 



This is how way we take into the account the inter-site correlations described in the proceeding section. 
Next we use these Green functions in the expressions for electron numbers 



n: 



00 



(38) 



where 



f(e) = 



1 + exp[£/fcBT] ' 
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Treating the self-energies as S^, which are solutions of Eq. (29), and using it in the first order approximation we can 
obtain 



2 

n 



= E P L W VL W = (U- F? n -zJ)- T o (±<r) A, (39) 



»=i 

where the antiferromagnetic energy gap A is given by the expression 

777 

A = F AF - (40) 
and Faf is the value of the effective exchange interaction given by 

F A F = Fin - ZJ, (41) 

with Fi n , the on-site exchange interaction (between different orbitals), being originally F® n given by Eq. (2). In the 
present case of a weak Coulomb correlation U (U/D <C 1), the on-site interaction, F; n , is augmented in the derivation 
process of Eq. (40), by the Coulomb interaction treated in the Hartree-Fock approximation. As a result, we have 

p. -f U + F ln [//£>« 1, 

ln ~\FP n for £//£>» 1. [A2 > 

It can be seen (Eq. (41)) that the inter-site exchange interaction, J, in an effective field (41) is opposing antiferro- 
magnetism. This is contrary to the case of ferromagnetism (see Ref. [10]). 

Using only the first order approximation for the self-energies (Eq. (39)) in Eq. (38) we obtain the following 
expressions for electron numbers: 



^ = »r = 2^ E \. p k a f^) + p rf(-E k )] , (43) 

where f(Ek) is the Fermi function 



2N 

k 



is the occupation probability of state (k, ±c), 

E k = y/el + A* (46) 

and 



The previous results for the antiferromagnetism (see Ref. [13]) are arrived at without the effect of bandwidth modifi- 
cation. In the present case, the kinetic energy still sums up to the previously used standard form (see Ref. [14]): 

+ tkS^) = ^> ( s * + s k a ) = 2 ]r s fc , ( 48 ) 



where now, in the expression for effective dispersion relation of the antiferromagnet; Ek = \/ z\ + A 2 , the dispersion 
relation s k was replaced by Sk = £°£af (n, m). 

Using definitions of the average electron numbers n^ a (7 = a,(i) (sec Eq. (22)), and expression (43) we obtain 
that in the mean-field approximation, the antiferromagnetic moment per atom (in Bohr's magnetons) is given by the 
following expression: 
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m = n a ~n a 



4 E p r) lf(E k ) f(-E k )]. (49) 



The chemical potential [i is determined from carrier concentration n on the base of the equation, which is also coming 
from the Eq. (43): 

k 

Inserting Eq. (45) to Eq. (49) and using relation (40) we obtain 

1 = -^af]T [/(£*) - /(-£*)]• (51) 

k k 



At the transition from AF to normal state condition (51) takes on the form 

_/ 

2e k 



1 = _^ f ^IM_M. (52) 



k 



The numerical analysis will be based on Eqs (52) and (50). It will show how the critical exchange field F^ F depends on 
electron concentration. The analysis of Eqs (51) and (50) will give the magnetization dependence on the temperature 
for different values of the inter-site interactions. 

5. Numerical examples 

To illustrate the numerical results for the antiferromagnetic ordering, the Eqs (51), (52) and (50) will be used in 
the integral form: 

1 = - f „/; ESLHi ,, m 



-D 

and 

D 



-D 



where variable e replaced the effective dispersion relation e k , E(e) — \J e 2 + A 2 and po is the constant DOS. This 
DOS for the weak correlation (U/D -C 1) is given by 

Po = { C > ° for other e (56) 
and for the strong correlation (U/D 3> 1) 



^D oL U A i 2 for - y/l - n/2D Q L AF < e < y/1 - n/2D L AF , 

P ° = ^ for - V^ d oLaf <e-U< ^/2D L AF , (57) 

for other s. 

The effective bandwidth D 7 which appears as the integral boundaries in Eqs (53)-(55), also depends on the Coulomb 
correlation U and for the weak correlation (U/D <C 1) is equal to 
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D = D L AF , (58) 



while for the strong correlation (U/D 3> 1) 



_ ( y/1 — h/2DqLaf for lower Hubbard sub — band, . . 

\ y / n/2D LAF for upper Hubbard sub — band. 

The parameter Iaf, which appears in the expression (19) for the bandwidth factor Laf, according to Eq. (5) depends 
also on the Coulomb's correlation U throughout quantity nf 7 . For the weak correlation (U/D <C 1) we have n\ a = 1 
in Eq. (5). Using additional values of from Eq. (22) we obtain 

Iaf = 2 "~" a ~ a f ■ (60) 
4(1 - m 2 ) 

For the strong correlation (U/D ^> 1), when nf 7 = 1 — (for the chemical level in the lower sub-band), we obtain 
from Eq. (5) that 

(2n - n 2 - m 2 ) (1 — n) 
(2 — ny — m 2 - 

when at m brings the following expression: 

rc(l-n) 

/o = ^7T- (62) 

This expression can be again calculated directly from the definition [4,9] of the average bond occupation with spin a 
in the paramagnetic limit of the band being split by Hubbard correlation 



OO 

J Po ^.f(e)de, (63) 



where the effective bandwidth D is given by (59) and po is given by (57). 

We assume the zero absolute temperature, for which the analytical expressions are obtained from Eqs (54) and (55) 
for the critical value of effective exchange field F^ F in the function of electron concentration n. In case of the weak 
correlation (U/D <C 1) , the effective exchange field F" F is given by 



per 

^AF — 



D 



jLaf- 



log 1 1 — n\ 

For the strong correlation (U/D ^> 1), the following expression is obtained: 

Do 



per 

^AF — 



2-n 



3n-2 



Laf- 



(64) 



(65) 



The critical value of the effective exchange field F A p depends strongly on the electron concentration and on the 
parameter of effective bandwidth Laf- Quantity F AF decreases with decreasing bandwidth. The components of the 
effective exchange field are the on-site and inter-site exchange interactions (see Eq. (41)). It is possible to notice that 
the on-site exchange interaction F ln favors antiferromagnetism, but the inter-site interaction J opposes it. 
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FIG. la FIG. lb 

FIG. 1. Dependence of the critical internal field on electron concentration for different values of the inter-site exchange 
interaction J in the case of weak (FIG. la) and strong (FIG. lb) Coulomb correlation, at to = 0.1 eV, z = 8, J' = J and 
V = At = 0; J = — solid line, J = O.lto — dashed line, J = 0.3to — dotted line, J = 0.5io — dotted-dashed line. 



Figs la and lb show the dependence of the critical internal field F-™ (given by Eq. (42)) on electron concentration 
for different values of the inter-site exchange interaction J. Additionally, the pair-hopping interaction J' = J [7] 
and V = At = is assumed. Analyzing the critical field for the weak Coulomb correlation (Fig. la) one can 
see that close to half-filling the increase in the inter-site exchange interaction J does not help antifcrromagnetism. 
The effect is reversed for smaller concentrations. There is a competition here between two effects. The increase of 
interaction J decreases the bandwidth Laf (see Eq. (19)) which decreases F^ F , but on the other side, it weakens the 
effective exchange interaction .Faf in accordance with Eq. (41). Similar behavior is observed for the strong Coulomb 
correlation (see Fig. lb) around the point of maximum in DOS which is shifted now due to the band splitting to the 
concentration of n = 2/3. 




FIG. 2. Dependence of the critical internal field on electron concentration for different values of the density-density 

interaction V in the case of weak (FIG. 2a) and strong (FIG. 2b) Coulomb correlation at to — 0.1 eV, z = 8 and 

J — J' = At = 0; V = — solid line, V = lto — dashed line, V = 3to — dotted line, V = 5to — dotted-dashed 
line. 
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Figs 2a and 2b present the dependence of the critical exchange field on the electron concentration for different values 
of the density-density interaction V. The other inter-site interactions are equal zero; J = J' = At = 0. Interaction 
V increases the bandwidth hence it increases the critical field F^ F . This in turn decreases the concentration range of 
possible antiferromagnetic state. 




FIG. 3a FIG. 3b 

FIG. 3. Dependence of the critical internal field Ff* on electron concentration for different values of the assisted hop- 
ping interaction At in the case of weak (FIG. 3a) and strong (FIG. 3b) Coulomb correlation at to = 0.1 eV, z = 8 and 
J = J' = V = 0; At = — solid line, At = O.lto — dashed line, At = 0.3t — dotted line, At = 0.5t — dotted-dashed line. 



The next interaction, which will be analyzed, is the assisted hopping interaction At. According to Eq. (19), the 
increase of At decreases the bandwidth, which decreases the total critical exchange field F£ F . Figs 3a and 3b show the 
dependence of the internal critical field JF^ on electron concentration for the weak and strong Coulomb correlation 
respectively, for different values of the assisted hopping interaction At. The other inter-site interactions are equal 
to zero; J = J' = V = 0. The change in the bandwidth Laf caused by assisted hopping interaction At, unlike 
the changes from all other inter-site interactions, is not symmetrical with respect to the half-filling. In effect, the 
dependence FF£(n) is not symmetrical with respect to the half-filling. Nonetheless, for U = oo we can draw curves in 
Fig. 3b only to the concentration of n = 1. The electronic states for n > 1 are not available for energetic reasons. 

1200 
1000 
800 
H 600 
400 
200 


0.8 0.85 0.9 0.95 1 

n 

FIG. 4 




FIG. 4. Dependence of the Neel's temperature (In) on electron concentration for different values of the inter-site interactions 
in the case of weak Coulomb correlation at to = 0.1 eV, z = 8 and Fin = 0.34 eV; J = J' = V = At = — solid line, 
J = J' = At = and V = 0.2t — dashed line, V = At = and J = J' = O.lto — dotted line, J = J' = V = and 
At = 0.15to — dotted-dashed line. 
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Fig 4 presents the dependence of Neel's temperature Tn, on electron concentration for different values of the inter- 
site interactions for the weak Coulomb correlation. All the curves are for the same internal exchange field, which is 
F; n = 0.34 eV. The curves with nonzero inter-site exchange interaction J and nonzero density-density interaction V, 
have smaller magnetization and the Neel's temperature. For half-filling (n = 1) they are 564 K and 936 K respectively. 
Without any inter-site interaction the Neel's temperature is higher, at half-filling it is In = 1000 K. The inter-site 
exchange interaction J decreases the critical temperature by decreasing the effective exchange interaction Faf (see Eq. 
(41)). This effect is stronger than the decrease of the bandwidth due to J, which increases Tn- The density-density 
interaction V decreases the critical temperature by increasing the bandwidth. Only the assisted hopping interaction 
At increases the Neel's temperature (to 1211 K at half-filling) by decreasing the bandwidth. 
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FIG. 5. Magnetization versus temperature for band filling n = 0.95, to — 0.1 eV, z = 8, Fi n = 0.34 eV and for various values 
of the inter-site interaction; J = J' = V = At = — solid line, J = J' = At — and V = 0.2to — dashed line, V = At = 
and J = J' = O.lto — dotted line, J = J' = V = and At = 0.15to — dotted-dashed line. 



Figure 5 shows the magnetization versus temperature for n = 0.95 and various values of the inter-site interaction. 
Again, all the curves are for the same internal exchange field F n = 0.34 eV. Similarly, as in Fig. 4 at constant value 
of internal field, the highest magnetization is obtained for nonzero assisted hopping interaction At. 
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FIG. 6. Magnetization versus temperature for band filling n = 0.95, to = 0.1 eV, z = 8, Fi n = 0.34 eV and for various values 
of the inter-site exchange interaction J at V — At — and J' — J; J = — solid line, J = O.lfo — dashed line, J = 0.12to — 
dotted line, J — 0.14to — dotted-dashed line. 
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Figure 6 shows the magnetization versus temperature for n = 0.95 and various values of the inter-site exchange 
interaction J. The value of internal exchange field is F; n = 0.34 eV. Magnetization is decreasing with increasing inter- 
site exchange interaction. For J w 0.122i the magnetization approaches zero as T — > 0, but for higher temperatures 
the reentrant antiferromagnetism occurs. For J > 0.15io the antiferromagnetic state at our values of n — 0.95 and 
F; n = 0.34 eV does not exist at any temperatures. 

6. Conclusions 

The most essential results of this paper will be recapitulated now. 

Introducing the inter-site interactions does not change qualitatively the results obtained much earlier (see [13]) for 
the antiferromagnetic ordering coming from the Coulomb correlation. In the case of the weak Coulomb correlation, 
there is a minimum of the critical exchange field at the half-filling. In the case of the strong Coulomb correlation, 
the minimum of these curves shifts close to the half-filling of the lower sub-band (n = 2 /3) and the half-filling of the 
upper sub-band (n = 4/3). 

The inter-site nearest-neighbor interactions lower the critical exchange field by the ratio of Laf = = 1 — 
rij^- + j- (—2 J — J' + V) Iaf- At the same time, the inter-site exchange interaction, J, increases the necessary critical 
exchange field by directly opposing it in the Weiss field; F" F = Ff* — zJ. The net effect of these interactions was 
illustrated in Figs 1-3. 

It has already been established that the superconducting cuprates have the superconducting band being split by 
the relatively strong Coulomb repulsion into two sub-bands, each with the maximum capacity of one electron. The 
new approach developed here allows for the antiferromagnetic instability in the strong correlation split band limit to 
occur at the half-filling, see Fig. 3b, when the assisted hopping interaction At w to, which reduces the bandwidth to 
zero, Laf ~ 0, at this concentration. At J = J' = V = we obtain from Eq. (65) the following formula: 



The other inter-site interactions (J, J', V) do not reduce the bandwidth at n = 1, since 1af(« = 1) = (see Eq. (61)) 
in the case of the strong Coulomb correlation. 

The reduction of the bandwidth to zero for strong correlation at half-filling in the purely itinerant model can also be 
caused by the limitation of electron jumps at this concentration, which according to Ref. [15] decreases the bandwidth 
t = toy/(l — n) I (1 — n/2). At close inspection of the expression for F AF coming from Eq. (54) at constant DOS and 
zero absolute temperature, the following formula is arrived at 



with Laf given by Eq. (19). This expression shows that F A p coming from this effect also tends to zero at the 
half-filling. 

In conclusion the antiferromagnetism at half-filling in the high correlation case can appear when the bandwidth 
goes to zero at this concentration. It can go to zero either as 1 — n(At/to), due to the presence of the assisted hopping 
(or some other interaction, whose strength is proportional to n) or as \f\ — n, due to the hopping limitation in the 
strong correlation limit. 

Both these situations can describe the high temperature superconducting cuprates, where at the half- filling and the 
strong Coulomb correlation there is, initially, the antiferromagnetic order. The superconductivity will appear upon 
doping only away from the half-filling at the electron concentration n ss 0.95 — 0.8 (see e.g. Ref. [16]). 




(66) 



£> Q £ A F\/(l-n)/(l-n/2) 
log(l — n) 



(67) 



13 



Appendix 



The stochastic potential in Hamiltonian (8) can be expressed as 

e = ( J + J')c+_ a c^ a + ( J - V)ct + At (n— + n~) . {AY) 

To arrive at this equation the stochastic value replaced each operator product in the square bracket of Hamiltonian 
(3). For example, the operator product; c- a C{ a was replaced by the stochastic value c7 a Ci a . It takes the values 1 or 
0. The probability of value 1, 7 CT , is the probability of electron with spin a hopping from the i to the j lattice site 
and back. The probability of the ii to j lattice site hopping is given by the product of probabilities that there is an 
electron with spin a on the i site and that the j site has empty states; nf (n* cr — n?)/ (nf + — nj) . The probability 
of forward and backward hopping is given by I a from Eq. (5) in the paper. 

According to this idea the stochastic potential i from Eq. (Al) takes on one of the following values with corre- 
sponding probabilities 



£l = 


J-V 






= / ff (l-/_ CT )(l-nr)(l-nT ff ), 


£2 = 


J + J' 




Pi 


= I^(l-Ia)(l-nr)(l-n-n, 


£3 = 


At 




pa 


= n- CT (l-/ (T )(l-/_ CT )(l-n7 CT ), 


£4 = 


At 




per 


= n7-(l-/ ff )(l-/_ CT )(l-nr), 


£5 = 


2 J + J' - V 




pa 


= 7 CT /_ ff (l-nT-)(l-nr), 


£6 = 


J-V + At 




Pi 


= 7 CT nr ff (l-7_ ff )(l-nT ff ), 


£7 = 


J-V + At 




Pf 


= 7 CT nT-(l-7_ ff )(l-nr), 


£8 = 


J+J' + At 




pa 


= I- a nr(l-Ia)(l-nJ°), 


£9 = 


J+J' + At 




pa 


= I- a nj°(l-I a )(l-n7 a ), 


£10 = 


= 2At 




pa 


= nj°n-°{l-I_ a ){l-I ry ), 


£11 = 


= 2J + J' - V + 


■At 


pa 


= l a I- a n~ a {\ - nJ CT ), 


£12 = 


= 2J+ J' - VI 


■At 


pa 


= I a I- a nJ°(l-n7 a ), 


£13 = 


-- J + J' + 2At 




pa 
^13 


= I- a n7 a n7 a {l-I a ), 


£14 = 


-- J - V + 2 At 




pa 


= I a n7 a n7 a {l-I- a ), 


£15 = 


--2J + J' -V + 2At 


pa 




£16 = 


-- 




pa 


= (l-7 (T )(l-7_ CT )(l-n- CT )(l-n7 CT ) 
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